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Classical diffusion of interacting particles 
in one dimension: General results and asymptotic laws 
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I consider the coupled one-dimensional diffusion of a cluster of A'^ classical particles with 

contact repulsion. General expressions are given for the probability distributions, allowing to obtain 

the transport coefficients. In the limit of large A'', and within a gaussian approximation, the diffusion 

constant is found to behave as A''"^ for the central particle and as (InA^)^^ for the edge ones. 

Absolute correlations between the edge particles increase as (InA'^)^. The asymptotic one-body 

distribution is obtained and discussed in relation of the statistics of extreme events. 
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I. INTRODUCTION 



' Low dimensionality, geometrical constraints and interactions between classically diffusing particles are expected 

to modify transport coefficients and/or the nature of the asymptotic regimes. As a trivial example, a brownian particle 



^ ■ subjected to a purely reflecting wall has a anomalous drift, its average position increasing as i^/^, (which corresponds 
to a vanishing velocity) whereas the centered second moment has a normal diffusive spreading with a lowered diffusion 



,; constant as compared to the diffusion without barrier. In the case of two particles with a contact repulsion, each of 
5t , them plays for the other the role of a fluctuating boundary condition which affects both the transport coefficients 
linked to the average position and the mean square deviation. The aim of this letter is to discuss such questions in 
the general case of N mutually interacting particles with a hard-core interaction. 

Classical diffusion with interactions does not seem to have drawn so much attention. A notable exception is the 
so-called "tracer problem" defined as the diffusion of a tagged particle in an infinite sea of other diffusing particles. 
The one-dimensional model was first solved by HarrisB and discussed in subsequent papersi'i. The main result is that 
the mean square dispersion of the position, Aa;^, displays a subdiffusive behaviour, which originates from the fact 
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that the motion of the tagged particle is, anywhere and at any time, hindered by all surrounding particles. In one 
dimension Harris found that Ax"^ grows as t^^^ at large times; this implies that the typical distance travelled by the 
tracer at time t goes like i^/^ instead of t^^'^ in the free case. More recently, Derrida et alfl found several exact results 
for the asymmetric simple exclusion procesJ^; recent bibliography and other results on such subject can be found in 
Mallick's thesiJl. 

The problem here considered is frankly different, although it belongs to the so-called single-file diffusion prob- 
lem encountered in many fields (one-dimensional hopping conductivityB, ion transport in biological mcmbranesi'i, 
channelling in zeolithcJl^) . At some initial time, a compact cluster of N pointlike particles is launched at the origin 
of a one-dimensional space; each of them undergoes ordinary brownian motion but has a contact repulsive interaction 
with its neighbours. As a consequence, the particles located at the egde of the cluster can move freely on one side and 
are subjected to a fluctuating boundary condition on the other, whereas the particles inside the cluster are subjected 
to such boundary conditions on either side. 

The questions to be solved are to find the (anomalous) drift of one particle and its diffusion constant as a 
function of the position within the cluster and of the number of N. In addition, two-particle correlations are worthy 
to be analyzed, as well as the asymptotic one-body probability distributions. 



II. ONE-PARTICLE TRANSPORT COEFFICIENTS 

At the initial time, the N particles are assumed to form a compact cluster located at the origin a; = 0, each of 
them having the same diffusion constant D as all the others. Due to the contact repulsion, two particles can never 
cross each other, so that order in space is preserved at any time; this means that the N coordinates Xi can be labelled 
so that: 

Xi < X2 < ■ ■ ■ < xn yt . (2.1) 
The solution of the diffusion equation for such an initial condition is the following: 

^ -xf/{4Dt) ^-1 

pixi,X2,...,XN;t) ^ mT\^^==^ T\Y{x,+i-X,) , (2.2) 

n=l nil 

where Y is the Heaviside unit step function {Y(x) = 1 if a; > 0, otherwise). In a recent and important paperlli, 
the general formal solution of the same problem with an arbitrary initial condition was given, using the reflection 
principle. ^From eq. (|2.2|), one readily gets the reduced one-particle density for the n**^ particle of the cluster: 



)1-JV 



(n - iy.{N -ny. 



1 + $ 



n-l 



1 - $ 



-x^/{4Dt) 



(2.3) 



In the last equation, $ denotes the probability integralEl, satisfying $(±oo) = ±1. The two factors (1 ± $) represent 
the steric effects on the n*'' particle due to the other ones. Knowing phK it is in principle possible to compute the 
first few moments giving the average position and the mean square dispersion for anyone of the particles. As a first 
result, one immediately observes that any moment of the coordinate has a quite simple variation in time, at any time, 



not only in the final stage of the motion. Indeed, using eq. (2^), it is readily seen that the the /c*'' moment < x^> 
increases at any time as t^^"^ since (DtY^^ is the only lengthscale of the problem. As a consequence, Vt, the average 
coordinate of the n}^ particle < Xn > (t) increases as t^^"^ - except for instance for the central particle of the cluster 
when N is odd -, whereas the mean square displacement Ax^ = < > — < Xn >^ increases as t. The drift, due 
to left-right symmetry breaking, is thus always anomalous and the diffusion is always normal. As a consequence, one 
can define, Vn and the following transport coefficients Vi/2.n ^nd 



<Xn>^ Vi/2,niN)t2 , 



(2.4) 



Axl = 2Dn{N)t , 



(2.5) 



It remains to find the functions Vi/2.„(A^) and Dn{N), which incorporate the dependence of the transport coefficients 
upon the mmrber of particles. For N — 2, one readily finds: 



< X2 > = - < Xi >= \/ -Dt , 



(2.6) 



and 



Axi = 2(1- -)t 



In the case of a non-fluctuating (fixed) perfectly reflecting barrier, one has: 



(2.7) 



4 2 
<x>=\l-Dt Ax^ = 2(1 )t 

" TT TT 



(2.8) 



Thus, for a cluster of two particles, each of them acting for the other as a fluctuating barrier, the drift is slowered 
and the diffusion is enhanced as compared to a fixed barrier. These facts are easily understood on physical grounds. 
Unfortunately, it does not seem possible to write the exact expressions of V^i/2, n I?„ in a closed, tractable 



form, starting from eq. (2_^). On the other hand, since it is worthy to analyze the case of a large number A'^ ^ 1 and 



3 



since the factors involving the $'s functions have rather sharp derivatives, especially when N is large, it is expected 
that a gaussian approximation can indeed produce the correct large- iV variation of the Vi/2,„ and Z)„. 

Let us first consider one of the two particles located at one extremity of the cluster, the right one for instance. 



From eq. (2.3), the one-body density can be written as (u = x/v4Di): 



1 



1 + 



TV 



(2.9) 



/iDt du 

This expression naturally has the form encountered in the statistics of extreme value Jll. When ^ 1, this is a very 
sharp function with a maximum Uq defined by: 

2^^ 



e "0 



N 



(2.10) 



Making now a gaussian approximation for p^^\x;t)^ one readily finds, up to logarithmic corrections: 



< xn > 



<Xi>'. 



N , 



2V^ 



(2.11) 



Thus, for large N , the coefficient for the anomalous drift has a logarithmic increase with respect to the number N of 
particles of the cluster: 



Vy2N{N) (X (ln7V)i/2 



(2.12) 



The fact that V1/2N increases with N is evident on physical grounds (all the "inside" particles are pushing on those 
which are at the edges), but this increase is extremely slow. In addition, the same approximation yields: 

e2/3 



N 



(27r)i/3ln ^ 



Dt 



(2.13) 



2v^ 



SO that: 



Di{N) = Dn{N) cx (InA^)-i . 



(2.14) 



Although the diffusion is normal, the diffusion constant decreases and tends toward for infinite N . In a pictorial 
way, the more there are particles pushing on its back, the less quickly spreads any one of the edge particles on either 



side of its average position, which drifts proportionnally to t^^"^. Note that, from eqs. ( 2.11 ) and ( 2.13 ), the relative 
fluctuations for the edge particles behave as (In A^)^^; this extremely slow decrease of fluctuations, as compared to 
iV~^/^ in ordinary cases, implies that convergence toward a large-number law, if any, is quite poor. 

The correctness of the gaussian approximation for the two flrst moments was checked by numerically computing 
the exact average position and exact mean square displacement and by looking at < xn > /[^Dt ln(A^/2^/7r)]^/^ and 



[Ax^/(4Dt)] In Fig. 1 displays the rather rapid convergence of such quantities toward constants at large N, 

confirming the validity of the gaussian approximation at least for the two first moments. 

Obviously, things go quite differently for the particle located at the center of the cluster, assuming N to be an 
odd number for simplicity. First, it does not move in the mean. Second, it must have a rather small diffusion constant 
as compared to the edge particles, since it is strongly inhibited by its numerous erratic partners on either side of it. 
Indeed, the gaussian approximation yields: 



Ax^jv+i)/2 ^ ^Dt . (2.15) 



This provides the large- A'' dependence of the diffusion constant for the central particle 

1 



D^n+dMN) (X - , (2.16) 



entailing that the fluctuations are now of the order of 1/viV. 

Thus, in any case, the diffusion is normal, as contrasted to the Harris' case for which Ax"^ oc t^^'^. Yet, note 
that in the N ^ oo limit, both Dn and -D(Ar-|-i)/2 vanish, which indicates a lowering of the dynamical exponent. The 
vanishing in all cases of the diffusion constants in the A^-inflnite limit signals the onset of a subdiffusive regime in the 
finite concentration situation. Considering the middle particle, which is surrounded by infinitely many other, this is 
in conformity with Harris' result. For the two edge particles, the marginal logarithmic decrease of Dn comes from 
the fact that the former still face a free semi-infinite space to wander in. 

Note that the scaling upon TV as described by eqs. (2.12) and (2.14) are the same as those obtained in ref.0; 



nevertheless, the asymptotic distribution law is not of the Gumbel type (see below). 

III. CORRELATIONS 

Statistical correlations inside the cluster are also worthy to analyze. As an example, let us consider the 
correlations between the two edge particles. For the latter, the two-body probability density is easily found from eq. 
( ^.2D as the following: 

p^^l{x^,XN;t) = ^^^^ MuN)-<i>{u,)f-'c-^-'^+<^Y{uN-u,) . (3.17) 

Two-body correlations are most simply measured by Ciat = < xiXn > — < xi >< xn >', making again a gaussian 
approximation for the two-body density, this correlator has the following approximate expression: 
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CiN{t) ^ 41n 



TV 



Dt 



(3.18) 



Due to scaling in space, the normalized ratio CiAr(t)/Ax^ is a constant in time. This constant turns out to be an 



increasing function of the number N of particles; from eqs. (2.11) and (2.13), one finds 



27r 



1/3 



In 



20F 



(3.19) 



Thus, increasing the number of inner particles enhances, although quite slowly, the correlations between the two edge 
particles. Far from inducing some kind of screening effect, repeated numerous collisions from inner particles enhance 
the statistical correlations between the edge particles. In a pictorial way, it can be said that the former act as "virtual 
bosons" by going from one to the other edge particles; the more they are, stronger is the effective (statistical) coupling. 

IV. ASYMPTOTIC DISTRIBUTION LAWS 



Interestingly enough, it is also posible to obtain the asymptotic form of the one-body distribution given by eq. 



(2.3). For the right particle {n — N), starting from eq. (2_^), one easily finds, with still u — x/\/ ADt > 



P^^(a;,t) 



N 



1 



1 
2^ 



N 



2U\/7T 



(4.20) 



V4:TrDt 

The maximum occurs for u ~ wq and the front is clearly asymmetric around uq - although the gaussian 
approximation, as shown above, well accounts for the large- dependence of the two first moments (expectation value 
and fluctuations). This asymmetry represents the pressure exerted by the inner particles on the edge ones. For the 
left particle, one simply has p^^-* (x, t) = p^^\—x,t). Fig. ^ shows that the large-iV expression, eq. ( 4.20| ), reproduces 



quite well the exact even for a moderately large value of N. From eq. (4.20), it is seen that p)^'(x, t) is not exactly 
a Gumbel distribution; on the other hand, the rescaled variable — ln(A^/2Y^) has, up to logarithmic corrections, 
the same dependence upon as a true Gumbel variable as far as the two first moments are concerned. 



As constrasted, starting again from eq. (2.3) for the central particle (n — {N + l)/2), the asymptotic form of 
the one-body density turns out to be simply the following normal law: 



1 



N 



-Nx^/(27rDt) 



TT V 2Dt 



(4.21) 



in agreement with eq.( ^.15| ). 
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Figure Captions 



1. Illustration of the asymptotic dependence upon N of the transport coefficients for the edge particles (average 
position, eq. (2.11) and mean square displacement, eq. ( ^.13 )). 



2. Comparison of the exact (solid line) and asymptotic distribution (dashed line) functions, respectively given by 



eqs. (p.3|) and (4.20), for a cluster of 1000 particles; the absissa is the reduced variable u ~ x/v'iDt. 
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